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(57) Abstract 



The invention concems a method for producing a Nyquist digital filter with null interference between symbols, designed to process 
a physical signal transmitted between a transmitter and a receiver via a transmission channel, said filter being a symmetrical filter P(z) = 
F2(z) of order N using an oversampling factor M=4, and forming an adapted pair comprising a transmitting filter and a receiving filter, the 
polyphase decomposition of F(2) being written: F(z) = Fo(z'*) + z-'FKz'^) + z-2F2(z4) + z-^FsCz^). The invention is characterised in that N 
is different from 4n, n being an integer, and: If N=4n+1, Fi(z)Fi(z) + z-»F2(z)F2(z) = 72""; If N«4+2, 2Fo(z)Fo(z) + F^i (z) + z-^F^^Cz) = 
7Z-"; If N=4n+3, Fo(z)Fo(z) + Fi(z)Fi(z) = 7Z-". F being the mirror symmetry of F and 7 being a non null constant. 



(57) Abr^g^ 

L' invention conceme un proc6d6 de realisation d'un filtre num^rique de Nyquist ^ interferences nulles entre symboles, destine ^ 
traiter un signal physique transmis entre un 6metteur et un rdcepteur via un canal de transmission, ledit filtre 6tant un filtre sym6trique 
P(z)=F2(z) d*ordre N mettant en oeuvrc un facteur de sur6chantillonnage M-4, et formant une paire adapt6e comprenant un filtre d*6mission 
et un filtre de reception, la decomposition polyphase de F(z) s'^crivant: F(z)=Fo(z'*)+z-'Fi(2'^)+2-2F2(z4)+z-3F3(z'*) caracteris6 en ce que 
N est different de 4n, n entier, et en ce que: Si N=4n+1, Fi(z)Fi(z)+z-*F2(z)F2(z)-7Z~"; Si N=4n+2, 2Fo(z)Fo(z)+F2i(z)+z-*F23(z)7r-"; Si 
N=4n+3, Fo(z)Fo(z)+Fi(z)Fi(z)=7Z-". F itant la symetrie miroir de F et 7 etant une constante non nulle. 
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Precede de realisation de filtres numeriques de Nyquist a 
interferences nulles entre symboles, et dispositif de filtrage 
correspondant. 

Le domaine de Tinvention est celui du filtrage de signaux numeriques, plus 
5 precisement rinvention conceme la realisation de filtres de Nyquist, tels que ceux 

utilises par exemple dans les systemes de transmission en bande de base et en 
monoporteuse. 

Les filtres de Nyquist jouent un role essentiel dans les systemes de 
transmission. Dans la plupart des cas, ce filtrage est reparti entre remission et la 
10 reception sous la forme de deux filtres dits en racine de Nyquist. A remission, le 

but est notamment de limiter la bande du signal emis. A la reception, le filtre doit 
eliminer les differents bruits pouvant alterer le signal utile, en particulier le bruit 
d' interferences dues a des canaux adjacents. Dans les systemes modemes, ces 
filtres sont realises sous forme numerique. II existe de nombreuses methodes pour 
15 leur synthese. Un element cle de cette synthese, outre bien sur le respect des 

specifications frequentielles, est Tobtention d'une interference entre symboles 
(lES) aussi reduite que possible. 

Pour assurer les meilleures conditions possibles de transmission, on 
recherche, en general, des filtres a coefficients reels qui, idealement, satisfont les 
20 criteres suivants : 

interference nulle entre symboles ; 
selectivite frequentielle ; 

filtres d'emisssion et de reception formant une paire adaptee ; 
linearite de phase. 

25 L' annexe 1 detaille ces differents aspects, et les contraintes 

correspondantes. Cette annexe, de meme que les suivantes, fait bien sur partie de 
la presente description a part entiere. 

Un systeme complet de transmission doit realiser generalement egalement 
d'autres fonctions indispensables telles que : 

30 - le filtrage analogique d' emission, qui a pour but de rejeter les 
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harmoniques en sortie du filtre numerique, ainsi que le filtrage 
analogique de reception, qui a pour fonction de limiter a la seule 
bande utile le spectre du signal a echantillonner ; 
les conversions numerique - analogique (CNA), a remission, et 
5 analogique - numerique (CAN) a la reception. 

L'invention conceme essentiellement Tobtention de filtres numeriques, et 
on ne detaille done pas les fonctions concemant la partie analogique. II convient 
cependant de noter que, en pratique, le filtrage analogique a mettre en oeuvre est 
nettement moins contraignant au sens de la raideur des filtres, que celui du 
10 numerique. Plus precisement, la complexite de realisation du filtre de rejection 

d'emision est en relation directe avec la valeur du facteur de surechantillonnage 
retenu. 

Quant aux convertisseurs, le fait le plus notable du point de vue du filtrage, 
est celui de la conversion numerique - analogique qui introduit un filtrage en 
15 sin(x)/x, qui pent etre compense classiquement par le filtrage analogique. Ce 

traitement pent egalement etre realise en numerique. 

Par ailleurs, pour des raisons de mise en oeuvre, les filtres sont 
generalement echantillonnes a frequence double ou quadruple. 

La figure 1 illustre la partie numerique d'un tel systeme de communication. 
20 Le signal a transmettre X(z) 1 1 alimente I'emetteur 12, qui delivre un signal filtre 

Y(z) 13 transmis via un canal de transmission 14 vers un recepteur 15 qui delivre 
le signal de sortie S(z) 16. Le signal X(z)ll est tout d'abord soumis a un sur- 
echantillonnage 121 puis a un filtrage 122 Fi(z). Dans le recepteur, le signal subit 
un filtrage de reception Fr(z) 151, puis une decimation 152. Avec T le retard 
25 elementaire lie au filtres numeriques d*emission et de reception, on a Ts = MT, ou 

Ts est la duree symbole et M le facteur de sur-echantillonnage. 

L'interet d'avoir un facteur de sur-echantillonnage eleve est notamment de 
faciliter la realisation du filtre analogique qui va suivre. Un inconvenient est que 
cela augmente d'autant la vitesse du convertisseur (CNA). En pratique, le 
30 compromis le plus raisonnable correspond a une valeur de M = 2 ou 4. 
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L'invention conceme specifiquement le cas ou le facteur de sur-echantillonnage est 
M=4. 

Des techniques de construction de filtres de Nyquist sont deja connues. 
Toutefois, les differentes families de filtres connues supposent toujours le 
relachement d'une des proprietes souhaitees parmi I'interference nuUe entre 
symboles, la linearite de phase (au moins a remission) et la paire de filtres 
adaptes. 

Souvent, dans les systemes de communication, on est conduit a accepter 
que le filtrage introduise une lES non nulle. Cela pose cependant parfois des 
problemes, par exemple dans le cas de modulations a tres grand nombre d'etats 
(64, 256,.,.) 

L'invention a notamment pour objectif de pallier ces inconvenients de Tetat 

deTart. 

Plus precisement, un objectif de l'invention est de foumir un procede de 
realisation de filtres de Nyquist a interferences nuUes entre symboles, et qui 
respecte egalement les conditions voulues de selectivite frequentielle et de linearite 
de phase, dans une configuration en paire adaptee. 

Un autre objectif de l'invention est de foumir un tel procede qui permette 
une mise en oeuvre pratique aisee, des filtres obtenus, et qui prenne en compte 
certaines contraintes pratiques, telles que les effets de la quantification des 
coefficients. Plus precisement, l'invention a pour objectif de fournir un procede 
garantissant le maintien de la propriete d'annulation de I'interference entre 
symboles apres cette quantification. 

Un autre objectif est de foumir un tel procede qui puisse etre mis en oeuvre 
a Taide d'une structure cascade, notamment sous la forme de treillis. 

Ces objectifs, ainsi que d'autres qui apparaitront par la suite, sont atteints 
selon Tinvention a Taide d'un procede de realisation d'un filtre numerique de 
Nyquist a interferences nulles entre symboles, destine a traiter un signal physique 
transmis entre un emetteur et un recepteur via un canal de transmission. Ce filtre 
est un filtre symetrique P(z) = F-(z) d'ordre N mettant en oeuvre un facteur de 
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surechantillonnage M=4, et formant une paire adaptee comprenant un filtre 
d' emission et un filtre de reception, dont la decomposition polyphase de F(z) peut 
s'ecrire : 

F(z) = Fo(z') + z-^F,(z')h- z-^-F,(z")+ z-'F3(z'). 

Selon rinvention, N est choisi different de 4n, n entier, et les coefficients 
de la decomposition polyphase de F(z) sont tels que : 

Si N=4n+1, Fi(z)f ,(z)+ z"'F2(z)f .(z)^ yz" 

Si N=4n+2, 2Fo(z)f o(z) + F;(z) + z'*F3(z) = yz"" 

Si N=4n-h3, Fo(z)f o(z) + Fi(z)P,(z) = Yz"" 

F etant la symetrie miroir de F et y etant une constante non nulle. 

Ainsi, le procede de Tinvention garantit, par construction, que TIES est 
parfaitement nulle. 

De fa9on preferentielle, N vaut 4n+3 ou 4n+l et : 

ledit filtre d 'emission effectue une interpolation par un facteur M=4 et 
presente une implantation correspondant a une decomposition polyphase 
dite de type II, telle que : 



F(z) = [z'^ z"' 



1] 



Fo(z") 



et ledit filtre de reception effectue une decimation par un facteur 
M=4 et presente une implantation correspondant a une 
decomposition polyphase dite de type I, telle que : 
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Les filtres ainsi obtenus sont simples a realiser et a mettre en oeuvre. 

Avantageusement, dans ledit filtre d'emission, on effectue une etape de 
filtrage suivie d'une etape d'interpolation par un facteur M=4. De meme, dans ledit 
filtre de reception, on effectue avantageusement une etape de decimation par un 
facteur M=4 suivie d'une etape de filtrage. 

Cette construction (structure permutee), permet de reduire d'un facteur 4 la 
cadence des operations. 

Selon un mode de realisation preferentiel de Tinvention, ledit filtre 
d' emission et/ou ledit filtre de reception presente une structure de realisation sous 
forme d'au moins un treillis. 

En effet, sous cette forme d 'implantation, la contrainte de reconstruction 
parfaite est integree structurellement. 

De fa9on avantageuse, ledit filtre d' emission et ledit filtre de reception sont 
constitues chacun d'une paire de composantes polyphases donnee respectivement 
par les equations suivantes : 




= gA(aJA(2)A(a„.,).„A(z)A(aJ q 



= gA(aJA(z)A(a„.,)...A(z)A(ao) 



1 a 



n 0 1 



avec : 



A(a) = 



et 



A(z) = 



0 z 



- 1 



-a 1 
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ou g est une constante non nulle de normalisation et a. des coefficients reels. 

Selon un mode de realisation parti culier de T invention, le procede met en 
oeuvre une structure a deux treillis. Selon une autre approche, il peut mettre en 
oeuvre une structure a treillis unique, fonctionnant a une frequence double. 

L'invention concerne egalement, bien sur, les dispositifs de filtrage 
obtenus a I'aide du procede decrit ci-dessus, II s'agit done de dispositifs de filtrage 
d'un signal numerique de Nyquist a interferences nuUes entre symboles, destine a 
trailer un signal physique transmis entre un emetteur et un recepteur via un canal 
de transmission, base sur un filtre symetrique P(z) = F^Cz) d'ordre N mettant en 
oeuvre un facteur de surechantillonnage M=4, et formant une paire adaptee 
comprenant un filtre d' emission et un filtre de reception, et dont la decomposition 
polyphase de F(z) s'ecrit : 

F(z) = F,(z^) + z-^F,(z")+ z-^"F2(2^)+ z-'F3(z"). 

Selon rinvention, N est different de 4n, n entier, et : 

Si N=4n+1, F,(z)f j(z)+ z'' F.Cz)? ^(z) = Y^'" 

Si N:=4n+2, 2Fo(z)Fo(z) + F;(z) +z"'F^(z) = yz"" 

Si N=4n+3, Fo(z)fo(z) + Fi(z)f ,(z) = Yz"" 

f etant la symetrie miroir de F et y etant une constante non nulle. 

D'autres caracteristiques et avantages de I'invention apparaitront plus 
clairement a la lecture de la description suivant un mode de realisation preferentiel 
donne a titre de simple exemple illustratif et non limitatif, et des dessins annexes 
parmi lesquels : 

la figure 1, deja commentee en preambule, illustre de fa9on 
schematique la partie numerique d'un systeme de communication ; 
la figure 2, commentee dans Tannexe 1, represente une 
specification frequentielle typique d'un filtre de transmission ; 
la figure 3, commentee en annexe 2, represente de fa^on generale 



un banc de filtres a 2 sous-bandes ; 

la figure 4 illustre un mode de realisation de la partie d'analyse 
d'un banc de filtres para-unitaires a 2 sous-bandes de la figure 3, 
sous forme de treillis ; 

la figure 5 rappelle la structure d'un systeme de filtrage emission- 
reception, dans le cas lineaire et para-unitaire, commentee par la 
suite ; 

les figures 6a et 6b representent, pour un filtre d'ordre 4n + 3, 
deux formes d'implantation du filtre d'emission de la figure 5, 
utilisant une decomposition polyphase respectivement en structure 
directe et en structure permutee ; 

les figures 7a et 7b illustrent, dans le cas de I'ordre 4n + 3, le filtre 
de reception de la figure 5, egalement dans le cadre de la 
decomposition polyphase, respectivement selon une structure 
directe et une structure permutee ; 

les figures 8a et 8b presentent le bloc treillis des filtres emission- 
reception des figures 6a, 6b, 7a et 7b, sous la forme 
respectivement d'un treillis direct et d'un treillis a sorties 
inversees ; 

la figure 9 illustre de fa9on detaillee un filtre racine de Nyquist 
d' emission d'ordre 4n + 3 a deux treillis, mettant en oeuvre les 
treillis inverses de la figure 8b ; 

la figure 10 presente un filtre de Nyquist de reception d'ordre 
4n + 3, mettant en oeuvre quatre treillis ; 

la figure lObis presente un filtre de reception mettant en oeuvre un 
treillis direct selon la figure 8a (voir annexe 4) ; 
la figure 1 1 rappelle les equations utilisees dans le cas du treillis 
inverse de la figure 9 ; 

la figure 12 presente un autre mode de realisation du filtre 
d'emission, mettant en oeuvre un treillis unique et f onctionnant en 
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cadence double ; 

la figure 13 presente un exemple de reponse en frequences 
optimisees d'un filtre racine de Nyquist d'ordre 43, selon 
rinvention ; 

5 - la figure 14 illustre la realisation d'un bloc Z, dans le cas d'un filtre 

d'ordre 4n + 2, tel que discute en annexe 3 ; 

la figure 15 presente le bloc d'implantation de la matrice P'^, 
egalement dans le cadre du mode de realisation de I'annexe 3 ; 
la figure 16 presente le schema complet d'implantation d'un bloc a 
10 4 entrees et 4 sorties pour Tordre 4n + 2, note Ma, et discute en 

annexe 3 ; 

la figure 17 presente un mode de realisation d'un filtre demi 
Nyquist d'emission pour I'ordre 4n + 2 ; 

les figures 18 et 19 (bloc N de la figure 18) illustrent un mode de 
15 realisation demi-Nyquist de reception pour I'ordre 4n + 2, discute 

en annexe 3 ; 

les figures 20 et 21 presentent respectivement les filtres racine de 
Nyquist d'emission et de reception pour I'ordre 4n+l (ou 4n4-5), 
tels que discutes en annexe 5 ; 
20 - les figures 22 et 23 illustrent les resultats obtenus a I'aide de la 

methode de synthese generique commentee en annexe 6. 
Comme deja mentionne, le procede de I'invention permet la realisation de 
filtres avec M = 4 verifiant I'ensemble des criteres suivants : 
lES nulle par construction ; 
25 - filtres a phase lineaire ; 

filtres emission et reception formant une paire adaptee. 
L' annexe 2 montre, apres avoir developpe 1' analogic entre une paire de 
Nyquist avec M = 4 et un banc de filtres orthogonaux a 2 sous -bandes, que deux 
situations peuvent se presenter, selon I'ordre du filtre. Plus precisement, on 
30 montre qu'il n'y a aucune solution pour les filtres d'ordre N = 4n, et que les cas 
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N = 4n +1 et N = 4n + 3 sont equivalents (th^oreme 2), en ce sens qu'ils 
produisent des composantes polyphases similaires. La condition pour obtenir une 
interference entie symboles nulle est pr6cisee par le theoreme 1. 

Le cas des filtres d'ordre 4n + 3 est d6crit par la suite, n s'avere Stre celui 
pennettant d'obtenir la realisation la plus simple. Le cas 4n + 2 est traite k V annexe 
3, et le cas 4n + 1 (ou 4n + 5) & Tannexe 5. 

D'une fa9on g6nerale, le schema de base mis en oeuvre est represent^ en 
figure 5. Pai rapport au schema initial g^n^ral de la figure 1, on note que Fr(z) 
122 et FR(z) 151 sont remplaces par F(z). II s*agit en effet d'une paire ad^t^e k 
phase lin^aire. 

On note egalement Tintroduction d*une iigne k retard z-r 51, oil r est une 
valeur entiere positive d^pendante de Tordre, qui permet dans tous les cas 
d' obtenir une IBS nolle avec un retard donn6. Le facteur multiplicatif g 52 pennet 
d* assurer qu*au retard pr^» la sortie soit bien identique k Tentrfe. 

Dans le cas N = 4n + 3, du fait des relations entre les differentcs 
composantes polyphases F(z) Texpression (20) de Tannexe 2 pent se r6&rire de la 
fagon suivante : 

F(z) = Fo(z*) + z-^/\(z') + z-'^a^)+z-'Fo(z*) (31) 

On sait qu'au delai pihs les filtres emission et reception sont identiques. 
Cependant, si Ton veut tiier le meilleur parti des traitements multicadence, c'est-a- 
dire interpolation a remission, d^imation k la reception, les schemas de 
realisation ne seront pas exactement les memes. 

A remission, ot Ton rdalise une interpolation par un facteur M = 4, les 
formes d'implantation representees aux figures 6a et fib, correspondent k celles, 
dites d&omposition polyphase de type H [4] (les references citees sont regroup6es 
dans r annexe 7), donn^e dans le cas de Tinvention par : 



FEUILLE DE REMPLACEMENT (REGLE 91) 

ISA/EP 
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F^(z^) 



F(z) = [z"' z 




(32) 



Dans Tune des representations (figure 6b) on notera la permutation des 
operations expansion (61) - filtrage (62), qui permet de realiser un traitement 
equivalent en travaillant a une cadence 4 fois moindre que dans le cas du schema 
initial (figure 6a), 

Dans la partie reception le filtrage est associe a une operation de 
decimation, pour laquelle une ecriture plus appropriee de F(z) est celle de la 

decomposition polyphase dite de type I [4], qui dans le cas present peut s'^crire : 



On liote que, par rapport h la figure 5, le d61ai est tel que r=L 

Deux variantes de Timplantation de cette demi^re equation sont donnees a 

aux figures 7a et 7b (symetriquement aux figures 6a et 6b). 

L' analogic entre le filtre F(z) d'ordre 4n + 3 et le banc de filtres d'ordre 

2n + 1, decrite par le theoreme 3 (annexe 2), permet aussi de re6crire Tequation 

sous la forme : 



1 




(33) 



z 



F (z) -F (z) 



(34) 



F (z) F (z) 



= gA(a )A(z)A(a )...A(z)A(a ) 

n n - 1 0 
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Les deux paires de coixfiosantes polyphases, miroirs Tune de Taucre, se 
deduisent done par les 2 equations qni suivent 



=gA(aJA(z)A(a„_,)...A(z)A(ao) 



(35) 



et 



=gA(a,)A(z)A(a„_, )...A(z)A(ao) 



(36) 



10 



15 



20 



Dans chacun de ces cas, la structure treiUis, representee aux figures 8a et 
8b, est la mdme. 

Les deux vaiiantes k sorties directes ou invers6es de la figure 8 pennectent 
ensuite de r6aliser les filties emission et reception suivant les schdmas lepr&ent^ 
rcspectivement aux figures 9 et 10. On notera que les symboles en sortie des filtres 
indiqurat seulement quels sont les filtres acti&. 

En utilisant la d&omposition polyphase du signal d*entree, on peut ensuite 
determiner un sch6ma equivalent a un seul treillis mais qui, dans ce cas, est 
contraint de fonctionner k la fiequence double. 

Dans le schdma propose on presente le treillis inverse sachant que, bien 
entendu, le treillis direct conduit k un schema de mSme type. La figure 1 1 rappelle 
les Equations de base de ce tceillis. 

Le signal de sortie du filtre d*6mission s*exprime par 

nz) = z-'[F^iz') + z-'F,{z') + z-'^(z') + z"'Fo(z*)]x(z") (37), 



ce qui peut se reecrire : 

Y(z) = z-po(z')-i-z-'F,iz*)]xiz')+z-'[F,iz')+ z-%iz')]K{z') 

= [r,(z^)+r'ii(z*)]2-' 



(38) 



FEUILLE DE REMPLACEMENT (REGLE 91) 

ISA/EP 
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Une division par 2 du facteur d'expansion 4 s ^ 2 1^ 2), puis une 
permutation de Fexpansion par 2 avec le filtrage produit des sorties filtrees telles 
que : 



On reconnait dans Yi(z) la sortie superieure du treillis inverse fonctionnant 
a cadence double, et, au signe du terme en facteur de z'^ pres. dans Y^(z) la sortie 
inferieure. Pour retrouver exactement la meme expression, il suffit done de 
multiplier cette sortie du treillis inverse par (-1)", 

Une realisation basee sur ce principe est represente a la figure 12. 
Pour la partie reception, on pent considerer deux cas : 

le cas general, avec 4 treillis, qui satisfont Tensemble des 
proprietes desirees (figure 10) ; 

le schema simplifie, pour lequel seul TIES nulle et la linearite de 
phase a remission sont garanties (figure lObis). 
Get aspect parti culier est di scute en annexe 4. 

Dans le cas de la figure 9, on fonctionne a la frequence la plus basse. Le 
cout de chaque filtre correspond a une structure treillis a n + 1 cellules, plus celui 
d'une inversion de signe. Pour obtenir un gain unitaire une multiplication par le 
facteur g est egalement a comptabiliser pour le schema du filtre de reception. 

Les calculs pour chaque cellule elementaire correspondent aux 2 
multiplieurs et additionneurs. Cependant, comme pour chaque treillis une des 
entrees est a zero, un multiplieur et un additionneur peuvent done etre supprimes. 
le cout est done de 1 multiplieur et 1 additionneur pour cette premiere cellule et du 
double pour les suivantes, c'est-a-dire (2n + 1) multiplieurs et (2n + 1) 
additionneurs. Pour chaque filtre nous aurons done (4n + 2) multiplieurs et 
(4n + 2) additionneurs. A la frequence la plus basse nous avons done a realiser 
(4n + 2) MPU et (4n +2) APU. 



Yo (z) = [ Fo (z' )X (z^- ) + z- * f , iz' )X iz' ) ] 
Y , (z) = [ F , (z^ )X (z^- ) + z- ^ Po )X (z^- ) ] 



(39). 
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Pour la partie reception, dans le cas du schema simplifie, illustre par la 
figure lObis, la complexite operatoire est d'un seul treillis plus celle des trois 
additionneurs, du multiplieur par g et de Finverseur. Les calculs effectues a la 
cadence la plus basse se chiffrent done a (2n + 3) MPU, (2n + 5) APU et une 
5 inversion. 

Dans le cas general illustre par la figure 10, a la frequence la plus basse, la 
complexite est de (8n + 4) MPU et APU, et une inversion. 

La methode de synthese se deroule par exemple en trois etapes, illustrees 
ci-dessous par le calcul d'un filtre d'ordre 43. 
10 - Etape 1 : A partir d'outils de synthese [3], il est possible d'obtenir un 

premier jeu de coefficients transversaux qui soit proche de 
TIES nuUe. Ce sont les coefficients initiaux de la table 1 ; 
Etape 2 : Par identification, nous pouvons calculer les coefficients 
treillis initiaux correspondants qui sont donnes sur la table 
15 2. Le jeu initial de coefficients transversaux n'etant pas 

exactement a lES nuUe, on obtient ainsi un nouveau jeu de 
coefficients transversaux, appele reconstruit (cf table 1) ; 
Etape 3 : Une optimisation locale des coefficients treillis initiaux 
permet en suite d'ameliorer un peu la reponse en frequence 
20 du filtre. Les tables 1 et 2 donnent les valeurs de ces 

coefficients treillis et transverses apres cette optimisation 
(tables 1 et2). 



25 



30 
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i 


initial 


Keconstruit 


Optimise 


0 


0.006302365568 


0.006302365568 


0.00374315552336009 


1 


-0.01335763186 


-0.01335763186 


-0.00764365688411938 


2 


-0.003586256877 


0.000139948824594686 


0.000256865997947587 


3 


0.0032S1590529 


0.000296616382849522 


0.00052452951560661 


4 


0.008254154585 


0.00825445268157875 


0.00210146572675043 


5 


0.008753151633 


0.00875251982827294 


0.00961384935212083 


6 


0.002954708412 


-0.000818038274240153 


-0.000565430770503378 


7 


-0.006948650815 


-0.00231664836261536 


-0.00210883853472157 


8 


-0.01474214438 


-0.0147620230725521 


-0.00994253905900314 


9 


-0.01369718742 


-0.0136538137768085 


-0.0139493173848354 


10 


-0.001682000235 


0.00199238622180159 


0.00140539009490696 


11 


0.015714602.54 


0.00939101589973353 


0.00785652590927169 


12 


0.026870884 


0.0269017853424573 


0.0234276140713132 


13 


0.021-12355083 


0.0209715158594175 


0.0216783222304717 


14 


-0.003340088297 


-0.0032097877874163 


-0.00339966101869637 


15 


-0.03511243314 


-0.0308000965511264 


-0.0274417538899884 


16 


-0.05.300379917 


-0.0530202725288858 


-0.0511658629696837 


17 


-0.0364661812 


-0.0360408764445379 


-0.0377795925012697 


18 


0.0221698768 


0.0140041435375816 


0.0161657062857809 


19 


0.1105730906 


0.116698398397514 


0.109490464868349 


20 


0.1994492114 


0.199658215352476 


0.199471282748299 


21 


0.2550483644 


0.253344105750997 


0.257329314844394 



Table 4: Coefficients transversaux des filtres symetriques de degre 43. 





Initiale 


Optimisee 




0.006302365568 


0.003743155523 


"1 


2.1194631945539 


2.04203561311232 




- 0.7583059650965 


- 0.71854479175763 




0.6811397085317 


0.73802350460893 


0C4 


- 0.4637577935126 


- 0.55360950120952 


ces 


0.4470434662649 


0.51415010432349 




- 3.7431288142180 


- 5.34790460147966 




- 1.5196264288798 


- 1.52189568702330 


as 


1.2422617481334 


1.44910973460487 




- 0.5612421724760 


- 0.63596605846521 


OlO 


0.1588040486787 


0.18869460160260 


an 


- 0.0222057611677 


- 0.06862284945003 



Table 2: Coefficients des structures treillis. 
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La reponse en frequence pour le filtre d'ordre 43 calcule de cette maniere 
est reportee a la figure 13. 

II s'agit ici d'un premier resultat, qui peut encore etre affine. Par ailleurs, 
bien sur, d'autres techniques de calcul peuvent etre envisagees, pour differents 
5 ordres et contraintes techniques. 

On rappelle notamment que le cas N = 4n + 2 est presente en annexe 3. Par 
ailleurs, 1' annexe 5 discute le cas N = 4n + 1 (ou 4n + 5). 

Des precisions sur revaluation de la complexite des filtres de reception et 
une methode de synthese generique sont developpes en annexe 6. On notera que 
10 cette methode de synthese s'applique aux trois types de solutions discutes ci- 

dessus. 
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ANNEXE 1 

La repartition du filtrage de Nyquist sur 2 filtres, de ce fait couramment 
appeles filtres en racine de Nyquist, ou en raccourci demi-N\^quist, doit tenir 
compte de contraintes liees a TIES, aux specifications frequentielles, a la linear- 
ite de la phase et a la propriete dite de paire adaptee. Dans ce paragraphe nous 
precisons ces different es notions en supposant, pour simplifier la notation, que 
r = 1. Par ailleurs, dans le cadre de la realisation proposee, nous ne traitons 
que du cas des filtres a Reponse Impulsionnelle Finie (RIF). 

1- Interference Entre Symboles (lES) 

Soit x(n) le signal d'entree et X{z) sa transformee en i.e. X{z) = 
J3„aj(n)2~", et convenons de noter suivant le meme principe Tensemble des 
signaux du systeme de transmission. La contrainte d'lES nuUe revient a im- 
poser que, en I'absence de bruit, le signal de sortie, note soit identique, 

a un retard pres, k X{z). 

Soit Y{z) le signal en sortie de Temetteur, il pent s'ecrire 

Yiz) = FTiz)X{z''). (1) 
En sortie du filtre de reception le signal, note U{z)^ a pour expression 
U{z) = FT{z)FR{z)Xiz''). (2) 



Le signal de sortie est done tout simplement 



1 ^^-1 
■"'^ fc=0 



^ FT{z^w'')FRiz^w'') 

k=0 



X{z), (3) 



OU w = est une racine Mieme de Tunite, La condition d'lES nuUe pent 
alors se resumer a 1 'equation 

± J2 FT{z^w^)Fn{z^w^) = z-^ (4) 

ou d correspond au retard introduit par les 2 filtres. 

Soit P{z) le filtre produit, i.e P{z) = Ft{z)Fr{z). Dans un formalisme 
habituel en traitement multicadence du signal, Tequation (4) pent alors se 
reecrire 

P{z) |m= z-''. (5) 
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Un filtre P{z) qui verifie cette relation est un filtre de Nyquist, ou encore un 
filtre Mieme de bande. 

Soient nr et ur les ordres respectifs des filtres RIF d 'emission et de recep- 
tion, le filtre P est done d'ordre np - nr + riR et i\ peut se mettre sous la 
forme 

P{z)=^f:p{n)z-\ (6) 

Si Tequation (5) n'est pas verifiee, on mesure habituellement la "distance^ 
par rapport a cette propriete par Tune ou I'autre des expressions qui suivent 

I PW I kMrtd 

OU 

P W kM^d 

2- Les specifications frequentielles 

La figure 2 represente une specification frequentielle typique pour les filtres 
en transmission numerique : 

• Le facteur de retombee, en anglais le "roll-off", definit pour ce filtres 
20 passe-bas, la bande passante par [O.ujp = ^{1 — p)] et la bande attenuee 

par [u;, = §{l + p),7r]. 

• Les specifications d'ondulation en bande passante, notee 5i, et en bande 
attenuee, notee sont en general les memes pour les filtres d^emission 
et de reception. 



15 



25 



30 



A la pulsation tt/.I/, on impose en general la contrainte Fji{7t/M) — 
FriTT/M) = x/2/2. 

3- La linearite de phase 

Le choix de filtres a caracteristique phase-frequence parfaitement lineaire 
est en general recommande pour les systemes de transmission numerique [1, p 
325]. Cela a en particulier Tavantage de preserver les instants de passage par 
zeros des trains binaires transmis. 
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4- La paire de filtres adaptee 

Dans le Ccls de modulations lineaires, si on suppose que TIES est nulle 
pour Tensemble du systeme de transmission, et que le bruit du canal est ad- 
ditif, blanc et gaussien (BBAG), on sait que la paire adaptee, i.e. Ft{z) = 
z~-^' Fr{z''^), avec N ^ iir = tir Tordre de chaque filtre, est optimale pour le 
critere signal a bruit (SNR) [2, pp. 51 sq.]. Le SNR qui vaut ^p^jj^jjrqp [3] 
atteint alors son maximum avec |IF/?l|^||Fr||^ = 1- 
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ANNEXE 2 

1- Bancs de filtres ortliogonaux a 2 sous-bandes 

Un banc de filtres a 2 sous-bandes peut se representer suivant le schema de 
la figure 3. De maniere a satisfaire la condition de non repliement les filtres 
de sj'nthese Go{z) et Gi{z) se deduisent directement des filtres d'analyse par 
Go(z) = Hii-z)etGiiz) = -Ho{-z). ' 

On peut montrer ensuite [4] que la condition de reconstruction parfaite 
(PR) peut s'exprimer a partir de la seule matrice polyphase 'R[Ho,Hi]{^) du 
banc d'analyse 



H[Ho,ifi](2) = 



H,,o{z) H^,r{z) 



(9) 



ou les HkA^) sont les composantes polyphases pour chaque filtre Hk{z)^ i.e. 
Hu{z) = j:z-^H,,{z^) (10) 

Dans le cas des bancs orthogonaux a 2 sous-bandes les filtres sont d'ordre 
impair 2n + 1 et les composantes poly-phases sont toutes de meme degre n. Le 
banc est alors PR si le determinant de la matrice polyphase est un monome, 
20 autrement dit si 



Det H(Ho,i/.](^) = HMH^^r{z) - H^^o{z)Hoa{z) = /3z'- , (U) 

avec (3 une constante non nuUe. 

25 Soit Q{z) le filtre miroir du filtre Q{z) d'ordre 7n, i.e. Q{z) = z'^^Qiz-'). 

Le banc orthogonal se caracterise par des filtres Hq{z) et Hi{z) dits conjugues 
en quadrature, ce qui peut s'ecrire Hi{z) = ±z'^^''"^^^Hq{z'''^). Dans ce qui 
suit, contrairement a la presentation faite en [4], nous choisissons Tecriture 
avec un signe plus, Dans ces conditions les composantes polyphases satisfont 

30 la relation 



(12) 
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Par consequent la matrice polyphase peut s'exprimer pax 



H[Ho.H,(0 = [ ^^ ^^^^ 

et la relation (11) peut se reecrire sous la forme 



(13) 



Det HiHoM^) = Ho,o{z)Hofl{z) + Hr,a{z)Hxdz) = 0z-" , (14) 

Le signal de sortie est alors tel que X{z) = 7r~(^'''+^)X(z), avec 7 une 
constante non nulla. La matrice 'H.[fj^,ff^-^{z) est paraunitaire et peut se mettre 
sous forme cascade avec [4] 



H[77o,H,](^) = pA(Q„)A(z)A(a„_i) • . . A(2)A(qo), 



(15) 



ou g est une constante de normalisation non nulle, qq, . . . , a'„ sont n+l nombres 
reels et ou les matrices A(a) pour a reel et A(r) sont definies par 



A(a) = 



1 Q 

-a 1 



, A(.-) = 



1 0 
0 



(16) 



Une realisation d'un banc de synthese base sur cette ecriture mathematique 
produit la representation treillis de la figure 4. 

2- Caracteristiques generales du cas de la paire adaptee a phase 
lineaire pour M = 4 

Pour la paire adaptee a phase lineaire nous avons forcement Ft{z) et Fr{z) 
sjaiietriques tels que Fri^) — Fn{z)^ que Ton peut noter avec la seule expres- 
sion F{z), Soit I'ordre de F[z), nous avons done 



N 



F{z) = Y:h^ 



-k 



(17) 



Le filtre produit P{z) s'exprime alors par 



2N 



P{z) = F\z)^Y.P(^)^-' 



Jb=0 



(18) 



avec, bien sur, p{k) = p{2N — k). P{z) satisfait la condition d'lES nulle (5) si 
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p{k) = 0 pour fc - AT = 4/ , Z / 0 , (19) 

On peut faire un premier constat en ce qui concerne les filtres d'ordre 
5 = 4n. 

Proposition 1 .- Soil F(z) un filire symitrique d'ordre An. F{z) ne peut 
tire a IBS nulle. 

Pour que F{z), filtre symetrique, soit a lES nulle il faut que fo i=- 0. Le 
monome de plus haut degre de P{z) est alors que son terme central est 

10 de degre 4n en 2"^ Comme la difference des degres est un multiple de 4, F{z) 

n'est pas a lES nulle. 

La decomposition polyphase de Fiz) pour M = 4 s'ecrit sous la forme 



15 



20 



25 



= Fo(2^) + z-^Fi(z^) + z'-'F^^z') + z-^F^{z^) . (20) 



a.vec 

ni 



n=0 



Pour realiser Tanalyse des differents cas possibles en fonction de Tordre A^ 
on utilise la propriete de symetrie de F{z) et le fait que 4n/ + / < 

• Si iV = 4n, le degre de Fo{z) est egal a n et celui des Fi{z),i = 1, 2, 3 est 
strictement inferieur a ii. D'apres la proposition 1 on sait que ce cas ne 
presente pas d'interet pour notre probleme. 

• Si N = 4n + 1, Fi{z) est de degre n, Fo{z) de degre < n tandis que ^2(2) 
et Fziz) sent de degre < n - 1. Fi{z) correspond alors par symetrie 
miroir a ^0(2). ^1(2) est done en fait de degre n. On a done 



30 



Fo{z) = F,{z) et F^{z) ^ F2{z) 



(22) 
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• Si N = 4n + 2, il vient 77.2 = n,7?i = n et 713 = 7^-1. Faiz) est le 
symetrique de ^0(^)5 on a done 

5 

F2{z) = Fo{z) , Friz) = F,{z) , Fs{z) = Fs{z) (23) 

• Si = 471 + 3, 4712 + / < 4n + 3, V/ (0 < / < 3), et les filtres Fi sont tous 
de degre < n. En effet, le terme de plus haut degre de Fo{z) est de degre 

10 4n + 3 donnant un terme de degre n a Fsiz). Le terme de plus haut 

degre de Fsiz) correspond pax symetrie au terme constant de Fo{z). On 
a done 

Fz{z) = Fo(^-) et F^iz) = A(^) (24) 

15 . 

Remarque.- Dans Tequatioii (22) on a ecrit a dessein Fo{z) = ^1(2) et 

non Fi{z) = Fo{z) car on est sur que le degre de Fi{z) est n puisque le terme 

constant de Fq{z) est non nul, mais pas que le degre de Foiz) est n. D'ailleurs 

ce cas, ou /i = /at-i = 0, est considere plus bas. 

On developpe alors le produit F^{z) en tenant compte de la decomposi- 
tion (20) et Ton obtient la caracterisation suivante des filtres F{z) verifiant la 
propriete d'lES nuUe. 

Theoreme 1 - Soit F{z) un filtre symetrique d'ordre N verifiant F{z = 
0) 7^ 0. Alors suivant les valeurs de N, F{z) est a IBS nulle si et seulement 

25 

• Si N = 47Z + 1, 

Fi{z)Fi{z) + z--'F2iz)F2iz) = 7-'-" , (25) 

30 • SiN = An + 2, 



2Fo(z)Fo{z) + F^z) + z-'Fl{z) = 7---" , 



(26) 
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• Si N = 4n + 3, 



10 



15 



Fo{z)Foiz) + Fr{z)Fy{z) = o'-^"" , (27) 

ou 7 est unt constante iion nulle. 

Le theoreme suivant niontre que toutes les solutions d'ordre = 4(n— 1)+3 
s'obtiennent a partir des solutions d'ordre A'' = 4n + 1. 

Theoreme 2 - Soit F{z) un filtre symetrique a IBS nulle d'ordre An + 
Ln > 1, alors si Fi{z),i = 0, ... ,3 designe $e$ composantes polyphases, Fq{z) 
est de degre strictement inferieur a n tt F^{z) pent s'ecrire sous la forme 
Fi{z) = z-^Ki{z). On a alors ki{z) = Fo{z) et le filtre F{z) de composantts 
polyphases [Fo{z), F2{z),F3(z),Ki{z)] est symetrique d'ordre 4{n - 1) -h S et a 
IBS nulle. 

Demonstration - Cela resulte immediatement de la relation (25) puisque 
le terme constant du membre de gauche est nul. On a done Fi{z = 0)Fa(z = 
0) = fofi = 0. Comme Fi{z) contribue avec fo au terme de plus haut degre, 
fo y^O et done fi = 0. En ecrivant Fi{z) sous la forme Fi(z) = z''^Ki{z) et 
sachant que Fi{z) — Ki{z), la relation (26), peut encore s'ecrire 



20 z-'K,{z)k,{z) + z~'F2iz)F2{z) = yz'^ , (28) 

Comme ki{z) = Fq(z) et done Ki{z) = Fo(-) on obtient la relation (27), 
pour un filtre ou Kiiz) se substitue a ^3(2:) et a done pour composantes 
polyphases : [Fo{z), F2{z), F3{z)Ju(z)], On peut verifier sur un exemple le 
passage d'un filtre d'ordre 9 a un filtre d'ordre 7. 

25 

Exemple 

F,{z) =f,+f,z''^+foz'^ =^z''{f, + foz-') =z-U<i{z) ^29) 

F2{Z) =/2+/3^-' 
F3{Z) =/3+M-^ 

30 Le filtre compose de [Fo{z), F2(z), F3{z)J<i{z)] est d'ordre 7 et s'exprime 

alors par 

fo + faz-' + f3Z-' + f4Z-' + Saz-' + hz-'^ + hz-^ + foz-' (30) 
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Un troisieme theoreme montre que les solutions d'ordre = 4n+3 s'obtiennent 
a partir des bancs de filtres orthogonaux a 2 sous-bandes. 

Theoreme 3 - Soit [Ho{z),Hi{z)] un banc de filtres orthogonaux a deux 
5 sous bandes et d'ordre 2n + 1 dont la matrice polyphase est donnee par (13), 

alors lefiltre H{z) dont les quatre composantes polyphases sont [Hq^o{z), Hi^q{z)j 
Hi,o{z),Ho^q{z)] est un filtre a phase liniaire, d'ordre 4n + 3 et a lES nulle. 
De fagon reciproque, tous les filtres de ce type s'obtiennent a partir des bancs 
de filtres orthogonaux a deux sous-bandes. 

Demonstration - La demonstration est immediate puisque les filtres a 
phase lineaire, d'ordre 4n + 3 et a lES nulle sont caracterises par la relation 
(27) qui, dans les conditions du theoreme devient equivalente a (14). 

15 



20 



25 



30 
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ANNEXE 3 
Paire de Nyquist d'ordre N = 4n + 2 

1- Factorisation initiale 

Le theoreme 1 en annexe 2 nous donne, entre autres, la condition que doit 
realiser un filtre F{z) d'ordre 4n + 2 pour etre a lES nuUe. On rappelle qu'elle 
s'ecrit 



2Fo{z)Fo{z) + F^{z) + z-'Fiiz) = 



(31) 



10 



15 



20 
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ou 7 est une constante non nuUe et ou les composantes F{{z) sont toutes de 
degre n, hormis F3{z), et satisfont les relations F2{z) = Fo{z), Fi{z) = Fi{z) 
et Fsiz) = hiz). 

Theoreme 4.- Soient a{,i = 1, . . . ,n + 1, n > 0, n + 1 constantes. On 
considere le vecteur colonne [Fo{z) , Fi{z) , F2{z) , Fz{z)Y defini par I'egaliU 



Fo{z) 
F,{z) 
F,{z) 
Fsiz) J 



= M(ai) Z M{a2) Z 



Z M(a„) Z 



1 

1 

0 



(32) 



Alors le filtre F{z), dont les composantes polyphases sont les Fi{z),i = 1, . . . , 4, 
est a phase lineairej monique, d^ordre An + 2 et a IBS nulle, 

Reciproquement, tout filtre monique d'ordre 4n + 2^ d phase lineaire et a 
IBS nulle admet une decomposition de cette forme. 

Demonstration.- Pour n = 0, les composantes polyphases sont les polynomes 
constants Fo{z) = 1, Fi{z) = F2{z) = 1 et Fz{z) = 0. Ces polynomes veri- 
fient trivialement la condition (31) ainsi que les conditions de symetrie et de 
degre. On a alors F{z) = 1 + aiz"^ + z~^ . Supposons maintenant que les 
composantes d'un filtre F{z) d'ordre 4n + 2, Fi{z),i = 0,...,3, verifient la 
condition (31) ainsi que les conditions de symetrie et de degre. Construisons 
Li(z)J = 0, ... ,3 par 



(33) 



r Lo{z) ■ 




r ^o(^-) ] 


L,{z) 


= M(ai) Z 


Fx{z) 
F7{z) 


1 L^{z) J 




[ F^iz) J 
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On a 



a? 



Loiz) = Fo{z) - ^z-'F^{z) + a,z-'F^{z) , (34) 

a? 
2 



9 

,2 



L,{z) = -a,Fo{z) - a,z-'F2{z) + (1 - ^)z-'Fz{z) , (35) 



L2(z) = -^Fo{z) + z-'F2{z) + arz-'Fs{z) , (36) 

L3{z) = {l + ^)F^{z). ■ (37) 

II est a peu pres immediat de verifier les conditions de degre et de symetrie 

sur les composantes Li{z),i = 0, ... ,3. En developpant Texpression analogue 

de (31) pour les composantes Li(z)^i = 0, . . . ,3, on obtient 

2Lo{z)Lo{z) + Ll{z) + z-'Ll{z) = 

il±^z'^ {2Fo{z)Fo{z) + F^{z) + z-'Fi{z)) ^''^ 

D'apres (31), la partie directe du theoreme est done demontree. 

La reciproque est demontree si Ton pent etablir que les relations (34-37) 
permettent de determiner de maniere unique un coefficient ai et les com- 
posantes Fi{z),i — 0,...,3 d'un filtre F{z) monique, a phase lineaire et a 
lES nulle. Comme on doit avoir Fo{z = 0) = 1 (puisque F{z) est monique et 
symetrique, son terme constant est egal a 1), d'apres (35) on a ai = — = 

0) . Si ai ~ 0, alors nous sommes dans le c£is particulier ou Lo{z) est de degre 
n, le terme de degre n+l etant nul, le polynome symetrique Li{z) a un terme 
constant nul et son terme de plus haut degre est nul, et enfin le terme constant 
de L2{z) est nul (mais pzts son terme de degre n+l dont le coefficient est 

1) . On inverse alors les relations (34-37) par Fo{z) == Lo{z), Fi{z) = L^iz), 
F2{z) = L2{z) quo z~^ et ^3(2) = Li{z) quo z~^ . On suppose done main- 
tenant que ai 7^ 0, mais ce n'est pas necessaire en fait. L'inversion formelle 
des formules (34-37) conduit a 

^'^'^ = (TfW {^"^"^ " "'^'^^^ " ^^'^^0 ' ^^^^ 

= TT^^3(^) , (40) 

= (2 ^a\Y (-f - ""'^'^^^ + ^''^ ' ^^^^ 

F;{z) = (2T^ («i^o(-') + (1 - ^f)L^{z) + a.L^iz)^ quo z'^ (.42) 
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II suffit alors de montrer que les quotients pax z"^ sont exacts. Comme F{z) 
est monique, Lo{z = 0) = 1, Li(z = 0) = -ai par definition de aj et, d'apres 
la relation (31), L2{z = 0) = -a?/2. Les membres de droites des egalites 
definissant F2{z) et F3{z) s'annulent done en ;^ = 0 et les quotients sont exacts. 
On pent alors verifier que les polynomes i = 0, . . . , 3 satisfont les conditions 
de sy metric adequates. 

2- Decomposition equivalente 

En notant par R(a') la matrice 



R(cv) 



cos a 
sin a 

0 

0 



— sin oe 0 0 

cos a 0 0 

0 1 0 

0 0 1 



(43) 



et en utilisant la formule (32), on obtient 

= R(ai) PZP^ R(a2) PZP^ 



PZP^ R(a„) PZ 



^o(^) 
L Fsiz) J 

Dans le schema de remission on note par Ma le bloc associe au vecteur de 
droite. On introduit alors la matrice Z definie par 



1 

1 
0 



Z = PZP^ = 



et compte tenu de 



1(1 + ^-') 
0 

|(-1 + ^ 
0 



-1 



) 



0 
0 
0 

1 



0 

1(1 + ^-^) 

0 



0 

0 
0 



(44) 



■ 1 ' 




■ v/2 - 








1 




0 


0 




0 . 



(45) 



on obtient, avec un facteur multiplicatif ^r, la decomposition 
Foiz) 

^i^S"^! = ^P^R(ai) Z R(a2) Z . . . R(a„) Z 
L F^iz) 



1 

0 
0 



(46) 



Complexite operatoire 
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L'equation (46) donne directement la structure de realisation des filtres 
demi-Nyquist d'emission et de reception. EUe correspond a un systeme a 4 en- 



trees et 4 sorties qui dans les schemas des figures 6 et 7 se substitue a Tensemble 
des 4 filtres polyphases. On notera que dans le cas du filtre d'emission les sor- 
ties sont a inverser par rapport a Tecriture (46) qui correspond, par contre, 
exactement, a celle du filtre de reception. Hormis cette diflference de detail, les 
structures pour chacun de ces filtres d'ordre 4n + 2 vont done correspondre a 

• n blocs matriciels de type Z et R(a/). 

• un bloc matriciel P"^. 

• un bloc vectoriel [1, a„+i, 0, 0]-^. 

• un multiplieur g, a ne comptabiliser que pour remission ou la reception. 

Examinons a present plus en detail le cout de realisation de chacun de ces 
elements. 

• Matrices R(a/), 

Les matrices R(a/),1 < I < n, sont des matrices blocs qui peuvent se 
mettre sous la forme 



ou I est la matrice unite et G/ la matrice de rotation pour Tangle ai, i.e. 



Le cout d'implantation est done tout simplement celui du treillis corre- 
spondant a la rotation G/, c'est-a-dire 4 multiplications et 2 additions. 




(47) 




(48) 



• La matrice Z. 



# 
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EUe permet d'associer a iin vecteur e d'entree, le vecteur s de sortie par 
la relation s — Ze, ce qui se traduit par les 4 equations qui suivent : 

^0 = \[{eo - 62) + ^~^(eo + 62)] 

52 = \[{e2-eo) + z-\eo + e2)] ^^^^ 

53 = es 

Le schema de realisation d'un bloc de ce type est represente a la figure 
14. On peut preciser que les multiplications par 1/2 peuvent se reduire 
a un simple decalage des donnees binaires. En plus de ces 2 decalages le 
cout est done simplement de 4 additions et d'une inversion de signe. 

La matrice P"^. 

EUe se presente sous la forme : 



/a 0 -a 0 \ 

0 10 0 

a 0 a 0 

V 0 0 0 1 y 



(50) 



avec a = Sa realisation peut egalement se deduire des equations 

d'entree-sortie, qui dans le cas present s'ecrivent : 

So = a(eo + 62) 

''^1 ^ (51) 

32 = a(-eo + 62) ^ ^ 

S3 = 63 

de maniere immediate on peut done deduire ensuite le schema de la figure 
15. 

Le cout operatoire est de 2 multiplications, 2 additions et une inversion 
de signe. 

Le bloc vectoriel [1, On+i, 0, 0]^, 

Son cout d'implantation est d'une seule multiplication. 
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Le schema complet d'implantation, conforme a I'equation (46), correspond 
done a un systeme a 4 entrees et 4 sorties, tel qu'on le represente a la figure 16. 
On retrouve naturellement ce systeme a remission et a la reception. L'ordre 
de complexite de chacun de ces filtres est done de 4n MPU et 6n APU. 

Ainsi que nous I'avons deja precise, a I'emission il est necessaire d'inverser 
les sorties. Dans le schema de realisation reporte a la figure 17, cette operation 
est symbolisee par le bloc J, correspondant a la matrice antidiagonale, i.e. 



J = 



/ 0 0 0 1 \ 

0 0 10 

0 10 0 

V 1 0 0 0 y 



(52) 



Le schema de reception pour l'ordre 4n + 2 est represente aux figures 18 et 
19. II se deduit des calculs qui suivent. 
On introduit la matrice Z definie par 



Z = 



i(l + ^-^) 0 \{l-z-') 

0 0 0 

1(1 -^-^) 0 \{\^z-') 

0 1 0 



0 

0 
0 



C'est la matrice obtenue en prenant les elements "miroirs" des elements de la 
matrice Z. On a 

Z Z = 2r , 

ou I4 est la matrice identite d'ordre 4. 
On verifie ensuite les identites 



25 



PP' =14 , R(a)-'=R(a)' = R(-a) 

pour toute valeur de a. 

En notant M le produit matriciel 

M = P^R(oi)Z...R(a„)Z 
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et N le produit matriciel 

N = ZR(-a„) . . . ZR(-ai)P 
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on a done 



NM = z-n. 



4 5 



10 



15 



et par consequent 



N 



r Foiz) 1 




- 1 ' 


Fi{z) 


= gz 


0 


L Fs{z) . 




0 . 



On introduit alors la matrice K definie par 



K = 



10 0 0 

0 10 0 

10 0 0 

0 0 0 1 



La figure 19 donne le schema de realisation du bloc correspondant a la 
matrice N et la figure 20 le schema de realisation du filtrage de reception qui 
en resulte. Si la constante g est egale a 1/(1 + a„+i), alors un signal transmis 
par le systeme d'emission de la figure 17 et regu ensuite par le systeme de la 
figure 20, produit un signal identique avec un retard de n + 2 echantillons. 
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ANNEXE 4 
Schema de reception a un seul treillis 

Si on suppose un fonctionnement dit "dos-a-dos", il est en efFet possible 
de simplifier le schema de la figure 10 en utilisant Tequation assurant une lES 
nulle, on obtient alors le schema a un seul treillis de la figure 10 bis. Dans ce 
cas seules TIES nulle et la linearite de phase du filtre d'emission sont assurees. 

La validite de ce schema de la figure 10 bis pent se verifier a partir de 
Tequation (37) du brevet, apres multiphcation par z~^. Pour un systeme dit 
dos-a-dos, le signal d'entree du filtre de reception est donne par 

Yiz) = z-'[Fo{z^) + z-'F^{z^) + + z'^ Fo{z')]X{z^) (53) 

En reference a la figure 10 bis, les signaux apres decimation sont, de haut 
en ba^, donnes par 



Y{z) U= z-'Fo{z)X{z), z-^Y{z) U= z-' F,{z)X{z), 



(54) 



z-'Y{z)U=^-'F^{z)Xiz), z-^Yiz) U= Foiz)X{z) ' 
En entree du treillis direct, les signaux sont done, de haut en bas, donnes 

par 

20 ^"M^o(^) + Fi{z)]X{z) et z'^iFoiz) - Fi{z)]Xiz). Le signal de sortie est 

done 

S{z) = 29z-'[Foiz)Foiz) + F^iz)F,{z)]Xiz) = 2gz-'[jz-% (55) 



avec 



7=n(i+^')- 

On verifie ainsi que pour p = I/7, le signal de sortie est identique a celui 
de Pentree au retard de n + 1 echantillon pres. 



30 
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ANNEXE 5 
Cas de Tordre = 4n + 1 

Remarque preliminaire : Pour conserver un treillis a n + 1 cellules, de 
coefficients avec 0 < i < n, tel que celui defini pour Tordre 4n + 3, nous 
present ons le cas de Pordre 4n + 5. 

Soient Li^i = 0,...,3 les composantes polyphases d'un filtre L{z) a lES 
nulle d'ordre 4n + 3. On peut lui associer un treillis construit comme dans 
la figure 8 a) du brevet avec un jeu de coefficients ai^i = 0,...,n. On a 
L2{z) = Lt{z) et Lsiz) = Lo(z). 

Le filtre F(z) dont les composantes Fi^i = 0, . . . , 3 sont definies par 

Fo{z) = Loiz), (56) 

F,{z) = z''Loiz), (57) 

F,{z) = L,{z), (58) 

F^{z) = Li(z), (59) 

est un filtre a lES nulle et a phase lineaire d'ordre 4n + 5. De plus le theoreme 2 
en annexe 2 du brevet etablit que de tels filtres sont necessairement obtenus de 
cette maniere. Bien entendu, la synthese d'un tel filtre est conduite de maniere 
a optimiser les coefficients a,- , i = 0, . . . , n pour obtenir de bonnes specifications 
frequentielles pour le filtre F(z), pas pour le filtre L{z) ! Cela conduit done 
a la realisation de schemas pour remission et la reception associee a F{z) a 
partir du treillis construit avec les coefficients ai et du treillis inverse associe. 

La decomposition polyphase de type I applique a ce filtre se traduit par 
Texpression 

25 F{z) = Lo{z') + z-'lz-^'Loiz')] + z-^Ly{z') + z-n,{z') (60) 

On peut verifier que, pour la partie emission, cette equation peut se traduire 
par le schema en treillis de la figure 20, dans laquelle on applique un retard 
donne par en sortie de F{z), 

Pour la partie reception, comme dans le cas de I'ordre 4n + 3, il est pos- 

30 

sible de realiser un systeme a 4 treillis, qui garantit Tensemble des proprietes 
recherchees. 
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De meme des simplifications apparaissent si on considere le systeme dit 
dos-a-dos. Le signal d'entree du filtre de reception pent alors s'ecrire 

Y{z) = [z-'Loiz') + z-'[z-'Lo{z')] + z-'Li{z') + z-«L,(^^)].Y(z^) (61) 

Pour recuperer, a un retard pres, le signal d'entree nous mettons alors en 
oeuvre le dispositif de la figure 21. 

En reference a la figure 21, les signaux apres decimation sont, de haut en 
bas, donnes par 

Yiz) 4-4= z-^Lo{z)X{z), z-'Y{z) U= z-' Lo{z)X{z), 
z-^Yiz) U= z-n^{z)Xiz), z-^Yiz) U= z-^Lriz)X(z) ' ^ 

De haut en bas, les signaux en entree du treillis sont respectivement donnes 

par 

z-^[Loiz) + L^iz)]X(z) et z'^[Loiz) - Li{z)]X{z). 

Nous trouvons done comme pour le cas de I'ordre 4n + 3, mais avec cette 
fois un retard de traitement de n + 2 echantillons, un signal de sortie qui pour 
g = 1/(27) ^st identique a celui d^entree. 
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ANNEXE 6 



1- L'evaluation de la complexite des filtres de reception 

Deux cas sont a considerer suivant que Ton veut ou non conserver Pensemble 
des proprietes recherchees. Ainsi pour les filtres d'ordre 4n + 3 et 4n + 5 

• Pour conserver TIES nuUe, la paire adaptee, la linearite de phase des 
filtres emission et reception, comme indique a la figure 10, 4 treillis sont 
utilises qui fonctionnent a la cadence la plus basse. Compte tenu que 
pour chaque treillis 2 multiplications et additions sont necessaires, sauf 
pour chaque premiere cellule (un multiplieur et un additionneur), la com- 
plexite operatoire est de (8n + 4) MPU et APU, plus une inversion de 
signe. 

• Pour un systeme dos-a-dos des simplifications sont possibles mais seules 
TIES nulle et la linearite de phase du filtre d'emission sont garanties. 
Pour les filtres de reception d'ordre 4?^ + 3 et 4n + 5 simplifies la com- 
plexite operatoire correspond alors a celle d'un seul treillis plus celle des 
3 additionneurs, du multiplieur par g et de I'inverseur. Les calculs efFec- 
tues a la cadence la plus basse se chiflfrent done a (2n + 3) MPU, (2n + 5) 
APU et une inversion. 

2- Une methode de synthese generique 

Par rapport a la methode de synthese proposee initialement nous avons 
actuellement une methode de synthese qui permet de satisfaire, pour les 3 
types de solution, une specification frequentielle telle que celle representee a la 



Pour un filtre F{z) de longueur donnee L, difFerente de 4n + 1, autrement 
dit pour un ordre diff'erent de 4n, et pour une valeur fixee du facteur p de 
retombee, nous cherchons a minimiser la fonction cout 



^{F) = snp{wp{l-\F{l)\)\wc[~-\F{e^m snp |F(0|2}. 



Nous utilisons Talgorithme FSQP ("Feasible Sequential Quadratic Pro- 
gramming'") developpe par I'equipe de A.L. Tits [5] pour la recherche du 



figure 2. 
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minimum d'un ensemble de fonctions cout regulieres non lineaires soumises 
a des contraintes generales, regulieres et non lineaires. Comme cette methode 
est une methode d'optimisation locale, le choix d'un point initial est neces- 
saire. Pour un poids donne w dans la bande d ' arret et un poids egal a 1 

5 a Porigine et en la frequence de Nyquist, c'est-a-dire 7r/8, nous calculous un 

filtre initial F*"^'* optimal pour la norme du minimax. Ce filtre F*^'* est en- 
suite utilise pour calculer directement un ensemble de coefficients treillis pour 
un filtre produisant une paire adaptee a lES nulle. Puisque F'"'* ne produit 
pas lui-meme une paire adaptee a lES nulle, I'identification entre les deux en- 

Q sembles de coefficients, les coefficients trans versaux de F^^^^ et les coefficients 

treillis, n'est pas exactement realisee. Nous obtenons par consequent un filtre 
different, note F''^*^ dont les coefficients treillis sont consideres comme point 
initial du probleme d'optimisation decrit par (63). Pour ce probleme, nous 
fixons wp = l^wc = 2,ws — 0.5. 

^ De plus, pour toute structure treillis optimisee pour une valeur donnee 

P05 nous utilisons egalement un algorithme de "continuation" pour obtenir 
une structure treillis optimisee pour une valeur p\ differente en considerant 
une suite de valeur intermedi aires pour p : le resultat de la structure treillis 
optimisee pour une valeur de p dans la suite est le point initial pour la valeur 
suivante de p. 

La figure 22 montre la meilleure attenuation dans la bande d'arret que nous 
avons pu obtenir pour une longueur et un facteur de retombee fixes. 

Considerons par exemple les valeurs L = 43 et /> = 0.5 : la meilleure 
attenuation, obtenue par optimisation directe, est egale a 52.28 dB. La reponse 
frequentielle du filtre d'emission correspondant, F^^*, est donnee par la figure 
23. Les coefficients transversaux et les coefficients de la structure treillis sont 
donnes dans la table 3. 
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Coefficients transversaux 



5 



10 



/o 


-1.033698 10-" 


/ii 


8.074275371845 10"^ 


/l 


-3.6954272174 lO"'' 


/l2 


1.5774357742194 10"^ 


/2 


6.60549905248 lO"'' 


/l3 


1.2603068512975 lO"'^ 


h 


-1.09979829792 lO"" 


/l4 


-6.358804372375 10"^ 


u 


4.31563793685 lO"" 


/l5 


-3.3366501346017 10"^ 


h 


6.08998879046 lO"" 


/l6 


-4.6538674170018 lO"-" 


fe 


6.39123817678 10"^ 


fl7 


-2.1044564354449 10"^ 


h 


-7.82626716253 10"'* 


fis 


5.1214636348857 10"^ 


fs 


-3.761395765875 10"^ 




1.51186813816588 10"^ 




-4.983816032565 10"^ 




2.39358058463701 10"^ 


fio 


-9.24070721312 lO"" 




2.74720061562686 10"' 



20 
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Coefficients treillis 




-3.1848708937 




2.648518499 10"' 


013 


-2.9385797132 


OLA 


-9.349457966 10"' 


«5 


9.246893837 10"^ 


0:6 


-1.308356726 10"' 


OLl 


3.6683584963 


^8 


-6.126257549 10"' 


^9 


-1.2465410017 




-1.439677875 10"' 




-3.5749582735 
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Table 3 : Coefficients treillis et tranversaux de I'exemple avec L = 43 at 

p = 0.5 (/42-. = /.). 
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REVENDICATIONS 

1. Precede de realisation d'un filtre numerique de Nyquist a interferences 
nuUes entre symboles, destine a trailer un signal physique transmis entre un 
emetteur et un recepteur via un canal de transmission, 

ledit filtre etant un filtre symetrique P(z) = F-(z) d'ordre N mettant en oeuvre un 
facteur de surechantillonnage M=4, et formant une paire adaptee comprenant un 
filtre d'emission (12) et un filtre de reception (15), 
la decomposition polyphase de F(z) s'ecrivant : 

F(z) = Fo(z") + z-^F,(z')+ z-'F,(z")+ z-'F3(z') 

caracterise en ce que N est different de 4n, n entier, 
et en ce que : 

Si N=4n+1, F,(z)f ,(z)+ z-^F,(z)f ,(z) = yz"" 

Si N=4n+2, 2Fo(z)fo(z) + Fi(z) + z-'F;(z) = yz'" 



15 



Si N=4nH-3, 



Fo(z)Po(z) + Fi(z)fiW = Yz" 



20 



f etant la symetrie miroir de F et y etant une constante non nuUe. 

2 . Procede selon la revendication 1, caracterise en ce que N vaut 4n+3 ou 

4n+l et en ce que : 

ledit un filtre d'emission (12) effectue une interpolation (121) par 
un facteur M=4 et presente une implantation correspondant a une 
decomposition polyphase dite de type II, telle que : 



F(z) = [5 



1] 



f,(z') 
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et ledit filtre de reception (15) effectue une decimation (152) par un 
facteur M=4 et presente une implantation correspondant a une 
decomposition polyphase dite de type I, telle que : 



1 




z 



-3 



3 . Procede selon Tune quelconque des revendications 1 et 2, caracterise en ce 
que, dans ledit filtre d'emission (12), on effectue une etape de filtrage suivie d'une 
etape d'interpolation par un facteur M=4. 

4 . Procede selon Tune quelconque des revendications 1 a 3, caracterise en ce 
que, dans ledit filtre de reception (15), on effectue une etape de decimation par un 
facteur M=4 suivie d'une etape de filtrage. 

5 • Procede selon Tune quelconque des revendications 1 a 4, caracterise en ce 
que ledit filtre d'emission (12) et/ou ledit filtre de reception (15) presente une 
structure de realisation sous forme d'au moins un treillis. 

6 . Procede selon la revendication 5, caracterise en ce que ledit filtre 
d'emission (12) et ledit filtre de reception (15) sont constitue chacun d'une paire 
de composantes polyphases donnee respectivement par les equations suivantes : 




roi 



fo J 



= gA(aJA(z)A(a„_,)...A(z)A(ao) 



1 



1 a 



[1 0 1 



avec : 



A(a) = 



et 



A(2) = 



0 z 



-1 



-a 1 
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oil g est une constante non nulle de normalisation et a. des coefficients reels. 

7 . Procede selon la revendication 6, caracterise en ce qu'il met en oeuvre une 
structure a deux treillis. 

8 • Procede selon la revendication 6, caracterise en ce qu'il met en oeuvre une 
5 structure a treillis unique, fonctionnant a une frequence double. 

9 . Dispositif de filtrage d'un signal numerique de Nyquist a interferences 
nulles entre symboles, destine a traiter un signal physique transmis entre un 
emetteur et un recepteur via un canal de transmission, 

base sur un filtre symetrique P(z) = F-(z) d'ordre N mettant en oeuvre un facteur 
10 de surechantillonnage M=4, et formant une paire adaptee comprenant un filtre 

d'emission (12) et un filtre de reception (15), 
la decomposition polyphase de F(z) s'ecrivant : 

F(z) = Fo(z^) + z-^F,(z")+ z--F,(z')+ z-'F3(z^) 

caracterise en ce que N est different de 4n, n entier, 
15 et en ce que : 

Si N=4n+1, Fi(z)f ,(z)+ z-'F2(z)f ,(z)=: yz"" 

Si N=4n+2, 2Fo(z)fo(z) + Fia) + z''F;(z)- yz'" 

Si N=4n+3, Fo(z)f o(z) + (z)f , (z) = yz"" 

f etant la symetrie miroir de F et y etant une constante non nulle. 

20 
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